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ABSTRACT 

The  change  of  variable 

t 

(1)  u(x,t)  =  /  0(x,T)dT 

0 

for  the  temperature  0  in  the  one-phase  Stefan  problem  leads  to  the  evolution 
inequality 

(2)  (ufc  -  Au  -  f)  (z  -  u)  >_  0,  for  all  z  >_  0  , 

where  u  ^  0  is  required.  This  inequality  is  to  hold  over  a  space-time 
domain  D  =  x  (0,T^)  with  a  Dirichlet  boundary  condition  imposed  on 

x  (0,Tq)  and  a  zero  initial  condition.  In  this  paper  we  examine  semi¬ 
discretizations  of  (2)  in  time  and  in  space  and  we  derive  the  respective  con¬ 
vergence  rates.  The  following  explicit  results  are  obtained: 

(3i)  || UM  -  u  ||  <  CAt  , 

IT  (D) 

where  l/1  is  the  H1 (ft) -valued,  piecewise  linear  in  time,  interpolant  obtained 
from  the  horizontal  line,  fully  implicit  Euler  scheme  applied  to  (2)  with 
At  =  Tq/M; 

(3ii)  ||U  -  u  ||  iCh2  , 

u  (D) 

where  U^  is  the  continuous  time,  finite  element  approximation  obtained  by 

applying  an  integrated  version  of  (2)  to  a  translate  of  the  finite  dimensional 
trial  space  of  C°  piecewise  linear  elements.  The  approximation  scheme  used 
to  define  U^  appears  to  be  new. 
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SIGNIFICANCE  AND  EXPLANATION 


Many  physical  phenomena  are  modelled  by  inequalities  rather  than  equations 


In  this  report  we  examine  a  dynamic,  or  parabolic  inequality,  which  character 


ice  melting,  under  certain 


which  separately  discretize  time  and  space.  The  rates,  as  well  as  one  of  the 


schemes,  appear  to  constitute  new  results  in  a  subject  which  is  still  rather 


undeveloped.  Subsequent  investigations  are  also  contemplated  for  inequalities 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


CONVERGENT  APPROXIMATIONS  IN  PARABOLIC  VARIATIONAL 
INEQUALITIES.  I:  ONE-PHASE  STEFAN  PROBLEMS 

Joseph  W.  Jerome 

§1.  Introduction. 

Parabolic  variational  inequalities  arise  in  various  modelling  problems  such  a? 
phase  transition  via  heat  conduction,  optimal  stopping  time  problems,  dynamic 
obstacle  problems  and  porous  medium  filtration  problems.  It  is  the  purpose  of  this 
series  of  papers  to  define  and  analyze  convergent  discretizations  for  the  afore¬ 
mentioned  problems  beginning  in  this  first  paper  with  the  one-phase  Stefan  problem, 
formally  equivalent  to  a  dynamic  obstacle  problem.  Although  it  might  be  preferred 
to  treat  all  of  these  problems  in  a  single  theory,  the  special  features  of  each 
make  the  unified  convergence  analysis  awkward  at  best.  In  order  to  avoid  subsequent 
repetition,  we  shall  restrict  ourselves  in  this  paper  to  the  analysis  of  semidiscre¬ 
tizations  in  space  and  time. 

The  one-phase  Stefan  problem  has  the  formulation  over  a  fixed  space-time 
domain  D, 

(i)  u  ^  0  , 

(1.1)  (ii)  u  -  Au  -  f  >_  0  , 

(iii)  (ut  -  Au  -  f)u  =  0  , 

subject  to  appropriate  initial  and  boundary  conditions  (cf.  Duvaut  [4],  Bensoussan 
and  Lions  [1],  Friedman  and  Kinderlehrer  [6]  )  •  The  specification  and  development  of 
this  formulation  will  be  given  later  in  the  introduction  together  with  a  statement 
of  the  physical  problem. 

An  exhaustive  numerical  analysis  of  (1.1)  and  other  variational  inequalities  has 
been  given  by  Glowinski,  Lions  and  Tremolieres  IS]  ,  though  rates  of  convergence,  and 
certainly  optimal  rates  of  convergence,  do  not  appear  to  have  been  derived  in  the 
literature.  For  (1.1),  we  give  a  horizontal  line  analysis  and  a  continuous  time, 
finite-element  analysis.  We  find  it  convenient,  however,  to  work  with  the  equivalent 
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formulation  (cf.  Theorem  2.4), 


(i)  u  >  0  , 

.2)  (ii)  (Tu)t  +  (u  -  r)  -  TAr  -  Tf  >_  0  , 

(iii)  t(Tu)t  +  (u  -  r)  -  TAr  -  Tf]u  =  0  , 

on  D,  where  r  specifies  the  boundary  values  of  u  and  T  =  (-A)  1  is  a  smooth- 

:  ig  operator.  The  semidiscretizations  are  defined  directly  in  terms  of  (1.2).  In 

‘  ms  of  the  corresponding  stationary  or  elliptic  formulation  associated  with  (1.1), 
this  amounts  to  minimizing 

(1.3)  E (u)  =  ^  (u,u)  -  (r,u)  .  -  (Ar  +  f,u) 

2  IT  (ft)  L  (ft)  H  (ft) 

2 

over  a  convex  set  in  L  (ft)  rather  than  the  usual  quadratic  energy  form  over  a 
convex  set  in  H^(ft).  Now  in  terms  of  the  horizontal  line,  discrete  time  approxima¬ 
tions,  the  use  of  (1.2)  merely  provides  an  efficient  way  of  estimating  the  rate  of 

2 

convergence,  0(At),  in  L  (D)  (cf.  Theorem  3.2).  In  this  case,  the  approximations 

based  on  (1.1)  and  (1.2)  are  identical.  However,  if  (1.2)  is  used  to  define  finite 

element,  continuous  time  approximations  taken  from  the  translate,  r  -  M^,  of  the 

C°  piecewise  linear  trial  space  with  homogeneous  boundary  values,  it  is 

2 

possible  that  the  L  projections  associated  with  (1.2)  are  not  identical  to  the 
H1  projections  defined  by  (1.1).  For  equations,  the  projections  coincide.  Thus, 
for  elliptic  equations,  replacement  of  the  energy  by  (1.3)  leaves  the  finite  element 
approximation  invariant,  so  long  as  the  discrete  H  *  norm  is  properly  defined. 

This  carries  over  to  parabolic  equations  (cf.  [10])  but  remains  an  open  question  for 
inequalities. 

In  fact,  the  previous  discussion  has  been  highly  simplified,  since  it  proceeds 
as  if  ft  were  a  polyhedral  domain,  compatible  with  the  vanishing  of  elements  of 

on  3ft.  Actually,  we  shall  assume  that  ft  is  smoothly  bounded  to  take  advantage 

of  the  known  regularity  of  u  in  this  case;  thus,  there  is  an  accompanying  boundary 

3/2  2 

layer  effect,  well-known  to  be  associated  with  h  convergence  order  in  L  .  We 


V'1 


'•  'i.Tft...  V,  i. 
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2 

shall  avoid  the  boundary  layer  by  using  the  L  formulation  of  (1.2)  to  define 
finite  element  approximations  indirectly  (cf.  (4.9))  in  terms  of  non-negative  definite, 
self-adjoint  approximations  of  T,  which  map  into  rather  than  M^,  i.e., 

the  trial  functions  are  unrestricted  on  3ft.  The  triangulation,  thus,  includes  non- 

2 

simplicial  elements  near  3.:.  Of  course,  this  removes  any  possibility  that  the  L 

projection  coincides  with  the  projection,  but  does  permit  the  derivation  of 

2  2 

0  (h  )  convergence  in  L  (D) ,  (cf.  Theorem  4.3),  provided  the  operators  T,  are 

h 

pointwise  non-negative,  which  we  explicitly  assume  (cf.  (4 . 8i ) ) .  This  always  holds 

for  T(cf.  Lemma  2.1)  and,  since  T.  is  defined  as  the  composition  of  T  and  the 

h 

H'*'  projection  E^,  always  holds  for  T^  in  one  space  dimension,  since  E^ 
coincides  with  the  interpolation  operator.  It  is  of  course  possible  that  the  point- 
wise  non-negativity  hypothesis  on  the  T^  masks  certain  boundary  layer  effects  in 
several  Euclidean  dimensions.  In  analyzing  the  finite  element  approximations  defined 
by  the  integration  method  just  outlined,  the  reader  should  not  confuse  this  method 
with  direct  methods  of  least  squares  type,  which  effectively  double  the  order. 

It  is  of  some  interest  to  compare  the  approach  to  the  one-phase  Stefan  problem 
using  variational  inequalities  with  that  using  the  classical  formulation 

(i)  e  =  &e,  e  >  o  , 

(1.4)  (ii)  0=0,  on  S  , 

30 

(iii)  A  cos  (v,l  )  =  |—],  across  S  , 

t  dV 

for  the  temperature  0  of  one  of  the  two  phases  in  a  two-phase  system  such  as  water- 
ice.  The  change  of  phase  results  from  heat  conduction  through  the  phase  ft  of 
positive  temperature,  which  is  separated  from  the  phase  ft\ft t  of  fixed  temperature 
0  by  the  moving  boundary;  S  is  the  time  profile  of  this  moving  phase  transition 

boundary,  A  is  the  latent  heat  of  fusion,  v  is  the  outward  normal  to 

30  30 

{(x,t)  :  0(x,t)  =  0}  and  [— 1  represents  the  discontinuity  of  —  across  S 

oV  dV 

directed  by  v.  The  fundamental  numerical  work  on  (1.4)  has  been  carried  out  by 
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Meyer  (cf.  [12])  who  devised  uniformly  convergent  free  boundary  approximation 
An  alternative  approach,  using  the  variational  inequality  formulation  (1.1), 
obtained  by  setting 


t 

u(x,t)  =  /  9(x,x)dT  , 

0 


sc':. cm' 


and 


f  (x) 


9  (x,0) 

-A 


X  £ 
X  i 


O 

"0 


ras  given  by  the  writer  in  [9] .  It  appears  that  the  two  approaches  are  equivalent 
with  respect  to  asymptotic  rates  of  convergence  for  the  computation  of 
S  :  {x  =  x(t),t}  in  one  space  dimension.  Indeed,  for  horizontal  line  semidiscre¬ 
tizations,  errors  in  computing  x  are  of  the  order  of  the  square  root  of  errors  ir. 
computing  u;  the  latter  is  computed  with  accuracy  At,  giving  a  net  accuracy  of 
i*St.  If  the  less  regular  variable  9  is  used,  the  error  in  computing  9  is  of 
order  /At,  yet  accuracy  in  x  is  proportional  to  accuracy  in  9;  the  net  result 
is  again  accuracy  of  order  /St.  This  heuristic  analysis  has  been  confirmed  in 
analytical  studies  of  the  one-phase  Stefan  problem  by  L.  Caffarelli. 

In  the  remainder  of  this  section  we  shall  discuss  the  precise  formulation  of  the 
variational  inequality  to  be  used  in  the  sequel.  It  can  be  shown  (cf.  [9])  that  (1.1) 
is  equivalent  to  a  pointwise  cone  inequality  (cf.  Abstract)  and  thence  to  the 
integrated  version  of  this  same  inequality.  More  precisely,  let  homogeneous  initial 
datum  and  a  boundary  datum  function  r  be  specified  satisfying 

(i)  r  e  X  =  W1'"(0,T0;L"(fi))  n  l“  (0,Tq;W2  (fi) )  , 


(1.5) 


(ii)  r(-,0)  =  0  ,  r  >  0  in  D  . 

Define  a  convex  set  C  by 

(1.6)  C={weX:  w  <_  r  in  D  ,  w (•  ,t)  c  H*  (R)  ,  0  <  t  <  TQ ) 

Finally,  let  f  be  given  satisfying 


(1.7) 


f  €  L  (D) 


Then  a  dynamical  obstacle  formulation  of  a  parabolic  variational  inequality  with 
solution  v  may  be  used  to  characterize  u.  It  may  be  written,  incorporating  initia; 
and  boundary  conditions,  as 
(i)  u  =  r  -  v  , 

(1.8)  (ii)  v  e  C  ,  v(-,0)  =  0  , 

(iii)  j  (v^  -  Av  -  fg)  (w  -  v)  >_  0,  for  all  w  e  C,  a.e.  in  (0,TQ) 

Here, 

(1.9)  fQ  =  r  -  Ar  +  f  . 


Remark  1.1.  The  formulation  (1.8)  is  convenient  for  piecewise  linear  finite  elements. 
It  is  immediately  seen  that  (1.8)  is  equivalent  to 

(i)  u  e  r  -  C  ,  u  ( •  ,0)  =  0  , 

(1.10) 

(ii)  /  (ut  -  Au  -  f)  (z  -  u)  >_  0,  for  all  z  e  r  -  C,  a.e.  in  (0,TQ) 
ft 

which  in  turn  is  equivalent  to  (1.1),  augmented  by  the  initial  and  boundary  conditions, 
and  to  the  non-integrated  form  of  (1.10). 

Remark  1.2.  Standard  methods  (cf.  Caffarelli  and  Friedman  [3])  show  that  (1.10)  has 
a  unique  solution  provided  f  is  Ca,  0  <  a  <  1,  in  each  of  connected  subdomains 
ftQ  and  ft^  with  common  smooth  boundary.  In  the  case  of  the  one-phase  Stefan  problem, 
ftp  and  ft^  are  the  regions  occupied  initially  by,  say,  water  and  ice.  In  [3]  it  is 
also  shown  that  8  =  u  is  continuous  on  D. 

Remark  1.3.  Throughout  the  paper  we  shall  assume  that  (1.10)  has  a  unique  solution 

u  e  r  -  C.  in  addition,  we  shall  require  some  further  regularity  properties  of  u  and 

ft.  The  hypothesis  on  u  is  described  in  (1.11)  below.  The  hypotheses  on  are 

basically  the  existence  of  classical  (Schauder)  regularity  theory  and  weak  maximum 

2 

principles  (cf.  Leninas  2.1,  2.2)  as  well  as  L  regularity  theory  (cf.  (2.2)).  We 
introduce  (1.11)  now.  Although  u  exists  only  as  a  measure  in  D,  it  is  to  be 
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x;  -acted  on  the  basis  of  explicit  modulus  of  continuity  estimates  in  [3],  coupled  with 


the  well-known  representation  (cf.  John  [111)  for  fundamental  solution  kernels 


es; cndina  to  -A,  that 


when  -A  is  viewed  as  an  operator  from  fCW)  onto 


The  idea  of  replacing  the  hypothesis  u  e  L  (D)  by  the  weaker  hypothesis 


was  introduced  by  the  author  and  M.  Rose  in  [10],  where  it  is  in  fact  proved 


(D)  for  the  discontinuous  enthalpy  H ( - ) 


>2.  Positivity  and  Equivalent  Formulation. 

Denote  by  T  the  inverse  of  -A  introduced  in  section  one .  Thus  T  is  an 
isomorphism  of  H  *  (ft)  onto  H*(r.).  Also, 


(2.1) 


<vv  -1 

H  (ft) 


-1  . 


<Tt1,t2> 


defines  an  inner  product  and  norm  on  H  (.2),  with  the  latter  equivalent  to  the 

standard  norm.  Here,  <  •,•>  denotes  the  duality  pairing  on  H*(ft)  x  H  *  (ft)  .  T  is 

2 

well-known  to  be  a  positive-definite,  self-adjoint  operator  in  L  ((2)  when  restricted 
to  this  space.  One  of  the  basic  hypotheses  on  ft  is  contained  in  our  assumption  that 
T  is  an  isomorphism  of  Hm  2  (ft)  onto  Hm(ft)  n  H*(ft)  for  m  ^  1: 


(2.2) 


|Tg| 


<  C||  g| 


m-2  . 


m  >  1 


H  (ft)  H  (ft) 

As  a  first  step  in  obtaining  the  formulation  (1.2)  we  state  the  following  familiar 
type  of  maximum  principle. 

Lemma  2.1.  T  has  the  positivity  property  in  the  sense  of  Korovkin: 

2  , 


(2.3)  g  e  L~ (ft) 

where  the  inequalities  hold  a.e.  in  ft. 


g  >_  0  •»  Tg  >_  0 


Proof:  Follows  directly  from  (1,  Theoreme  5.1,  p.  83]  or,  alternatively,  from  the 

00  »  2 

classical  maximum  principle  [7,  p.  35]  applied  to  the  dense  subset  C  (ft)  of  L  (ft). 
Lemma  2.2.  Suppose  u  satisfies  (1.10).  Then  u  satisfies  the  initial  condition, 

1  (2.4)  u ( • ,0)  =  0  , 

and  for  all  0  <  t  <  Tq,  the  boundary  condition, 


(2.5) 


and  the  integrated  inequality 

(i) 


u(  • ,  t)  =  r  ( •  ,t)  ,  on  3ft  , 


u  >  0,  in  ft 


(2.6) 


(ii)  (Tu>t  +  u  -  TAr  -  r  -  Tf  >_  0,  in  ft  , 

(iii)  /  [ (Tu)  +  u  -  TAr  -  r  -  Tf]u  =  0  . 

ft  C 
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Proof:  (2.4),  (2.5)  and  (2.6i)  are  immediate  and  (2.6ii)  follows  fro-  •  —  . 

by  applying  T  to 

(2.7)  u  -  Au  -  f  >  0  . 

Note  that  here  we  use  the  equivalent  characterization  (1.1).  Now  note  that  by 
(2.6i,ii)  it  suffices  to  prove 

(2.8)  /  [(Tu>  +  u  -  TAr  -  r  -  Tf]u  <_  0  . 

SI 

Thus,  if 

(2.9)  S2^={x£j2:u(x,t)>0}, 

then  12^  is  open.  Suppose  t  is  fixed  and  (1*  is  any  ball  in  .  ^  satisfying 

S2.  c  (2*".  In  particular,  u(x,t)  >  c  >  0  for  x  e  12..  For  almost  all  t  c  (0,T.)  , 

*  +  -  * 

1  ~  2+a  - 

u(-,t)  e  C  (S2)  so  that,  by  the  Schauder  theory  [7],  Tu(-,t)  e  C  (".)  ,  some 

0  <  a  <  1.  In  particular,  if  <p  e  Cq(S24),  i(/  ^  0,  there  is  an  e  >  0  such  that 

eiJiT (u ( •  ,t) )  <  c,  on  (2^ 

Thus,  the  choice 

z  =  u  -  eiJiTu 

leads  to  z  e  r  -  C  and 

/  (u  -  Au  -  f)\l>Tu  <_  0 

n  t 

2 

Since  ty  is  arbitrary,  we  may  set  ij;  =  ill  where  <Ji  -»■  1  in  L  (  '.)  .  This  leads  to 

n  n 

/  (u  -  Au  -  f)Tu  <  0 
(2  * 

and  then  immediately  to  (2.8)  for  almost  all  t  e  (0,TQ) ,  and  hence  for  all  t  by 
standard  continuity .  • 

Lemma  2.3.  Suppose  u  satisfies  (1.10).  Then  u  satisfies  (2.4),  (2.5)  and  (1.2) 

for  all  0  <  t  <  T  . 

0 

Proof:  It  is  enough  to  prove  that 

(2.10)  (T(ut  -  Au  -  f ]  }u  <_  0  in  52^  ,  0  <  t  <  TQ  , 
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where 


is  given  by  (2.9)  . 


It  follows  from  (2.6)  that 


(2.11) 

for  all 


/  (T [u  -  Au  -  f)  }(z  -  u)  0 

a 


Z  =  U  +  111 


supp  wca  ,  0  <  Z  €  L  (.':) 


If  is  any  compact  subset  of  SI  and  z  =  u  -  euy  .  0  <  e  <  1, 

+  * 

(2.11) , 


0  ^  E  /  {T[u  -  Au  -  f) }u\0 
S2  1 

so  that  (2.10)  holds.  ■ 

Theorem  2.4.  The  solution  of  (1.2),  (l.lOi)  is  unique.  In  particular 
the  solution  of  (1.10). 

Proof:  If  u1  and  u2  are  solutions  of  (1.2),  (l.lOi),  then 


3u 

{T[  IT  “  4ui  ~  f!}(ux  -  u2>  i  0 


3u2 

fT I  -  Au2  -  f)  }  (Ul  -  u2)  >_  0  . 


Subtraction  and  integration  over  ft  gives 


9(U1~U2) 

(  '  ux  ■  u2)  _x  +  <U.  -  u  ,  u  -  u  )  ? 

3t  1  2  H  l(0)  12  12  l2(!2) 


<  0 


Integration  in  t  gives 


(2.12)  |  llux  <  - -t)  -  “,(•  -t)  ||  2  +  /  ||  u  -  u  ||  2  ds  < 

H  ^  (SI)  0  1  2  l/(S)) 

for  0  <  t  <  T.  Uniqueness  is  immediate  from  (2.12)  and  the  remaining 
follows  from  Lemma  2.3. 


we  have  fror 


it  is  given  b 


0  . 


statement 


Remark  2.1.  Lemmas  2.2  and  2.3  form  the  basis  of  a  circle  of  equivalences  for  the 
solution  u  of  (1.10)  as  in  section  one.  We  note  for  future  reference  that  the 
solution  of  (1.10)  may  be  characterized  by 


C.13) 


(i)  u  t  r  -  C 


u  (•  ,0)  =  0 


(ii)  /  [  (Tu ) t  +  (u  -  r)  -  Tir 


for  all  0  <  t  <  Tq  and  all  z  £  (  )  ,  z  0. 


i 


I 


i 


1 
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-  Tf]  (z  -  u)  s 


;  3 .  Horizonal  Line  Analysis. 


Definition  3.1.  Let  M  be  a  positive  integer,  let  it  =  T./M  and  set  t  =  rait, 
-  0  m 

0  <  m  <  M.  With  u„  =  0  the  given  initial  datum,  let  {u  }M  „  be  the  recursively 
—  -  0  ’  m  m=0 

generated  sequence  of  solutions  of  the  elliptic  variational  inequalities 

(i ) 


u  £  r  -  C 
m  m  m 


(3.1) 


u  -  u  , 

(ii)  /  (  -2 — -—2 - Au  -  f  )  (*  -  u  )  >  0,  for  all  •f  £  C 

'  At  m  m  m  —  m 


obtained  formally  from  (1.7)  by  the  identifications  u(-,t  )  ~  u  and 


—  (t  )  ~ 

dt  v m' 


m-1 


At 


Here 


(3.2) 


(i)  r  =  r(-,t  ) 
m  m 


(ii)  C  =  (U'  £  W2'  (ft)  :  tji  <  r  } 
m  —  m 


f  =  f(*,t  ) 
m  m 


Remark  3.1.  A  unique  solution  u  e  W  (ft)  of  (3.1)  is  k^own  to  exist,  provided 
■  m 

f ( * , t )  is  piecewise  Ca  in  subdomains  of  ft  independent  of  t  and  separated  by 
a  smooth  boundary,  as  described  in  the  introduction.  Although  the  result  of 
Frehse  [5]  assumes  a  global  Ca  property  for  f(-,t),  modifications  (cf.  (9)) 
show  that  the  piecewise  result  is  true.  This  is,  of  course,  the  case  of  physical 
interest.  Appropriate  modifications  of  the  result  of  Brezis  and  Kenderlehrer  [2] 
also  yield  the  result.  The  technical  difficulty  involves  the  global  boundedness  of 
the  second  derivatives. 

Lenina  3,1.  The  relation 
Tu 


(3.3) 


Tu  i 

/  t  — m  .  "I~1  +  u  -  TAr  ( •  ,t  )  -  r  -  f  ]u 

'  At  m  m  m  m  m 


holds  for  all  1  m  M. 

Proof:  This  proceeds  by  a  straightforward  repetition  of  the  proof  of  Lernna  2.2. 

Theorem  3.2.  There  exists  a  constant  C,  independent  of  At  and  given  by  (3.13) 


below,  such  that 


•  r. 


(3.4)  max  ||u(-,t  >  -  u  ||  _  +  J. 

0<m<M  H  (n)  m=0 


+  l  ll“(-'V  -  Um  ll2f  At  1  cut)'  . 

m=0  L  (;■) 


In  particular,  if  u  is  the  step  function 

s ,  it 


u  (x,t)  =  u  (x)  ,  x  £  !)  ,  mit  <  t  <  (m  +  l)it  , 

s  /  At  m  — 


for  0  <  m  <  M  -  1,  then  the  estimate 


;!u  -  us  At  II  2  i  c(At) 

S'fl  L  (D) 


holds  for  some  constant  C. 


Proof:  Setting  t  =  t  in  (2.6iii)  gives,  since  u  >0, 
m  m  — 

(3.7)  /  l  <-.t  )  +  u(- ,t  )  -  Tir(-,t  )  -  r  -  f  )  (u(-,t  )  -  u  )  <  0 

^  at  m  m  m  mm  m  in 

Similarly,  we  obtain  from  (3.3)  the  reverse  inequality, 

Tu  -  Tu 

(3.8)  /  (  ■  ”  - -r ~  +  u  -  Tir(-,t  )  -  r  -  f  1  (u(-,t  )  -  u)  >  0  . 

*  At  m  m  mm  m  m  — 


For  simplicity,  we  set 


ym  =  u(-,t  )  -  u_  • 
m  in  in 


Then,  subtraction  of  (3.8)  from  (3.7)  yields,  after  an  appropriate  addition  and 


subtraction. 


Tu(-,t  )  -  Tu(.,t  > 

i  /t  - ^  V’ym 


r,  J<r  'Tw-t'V!^ 


■  f&r  (:.a)dadT}yn 

(2  t  ,  t  3t 

m-i 


m  mL2rn 


where  z  is  defined  by  this  equation.  Multiplying  through  by  .‘t,  making  use  of 
m 


the  inequality 
(3.10) 


I  <w-z>  I  <  \  llwll 2  +  \  11*11 2 


and  summing  on  m  =  1,  .  .  .  ,  k  yields 

k  .  .  k 

l 


(3 


i  m;.i. 


1*1  ?  At 
m  L2  C  ) 


H  ((2)  m=l  L  (ft)  m=l 

It  remains  to  estinate  the  right  hand  side  of  (3-11).  We  have,  by  the  Cauchy-Schwarz 
inequality  and  the  limiting  integral  form  of  the  triangle  inequality. 


(3.12) 


l’mll2±'ir/tm  r  11^  (-.0)11%  d,dT,2 

”  t  ,  T  3t  L  (SI) 

m-1 


<  At 


32Tu,,  2 


2  I'  2  2 

3t  L  ft  .  ,t  ;lT(f2)) 
m-1  m 


so  that  (3.4)  follows  with 
(3.13)  C 

Finally,  (3.6)  follows  via  the  triangle  inequality  from  (3.4)  and  the  estimate 


32Tuii  2 

2  II  2 

3t  h  (D) 


(3.14) 


M-1  tm+l 

_  M-1  ^+1  t 

I  /  /  |u<-.t)  - 

u( -#tm) 

I2  =  I  I  III 

m=0  t  0 

m 

m=0  t  fit 

m  m 

M-1 

Sn+l 

i  At  l 
m=0 

I 

llutl|22  2 

L  (t  ,t  , ;L  (ft)) 
m  m+i 

3u 

3t 


=  At2  Hu  "2 


t  "  2 
l/(D) 


The  theorem  is  now  proved 


:: 4 .  Continuous  Time  Finite  Element  Approximations. 


For  h  >  0,  let  be  a  triangulation  of  the  given  domain  ft.  Thus, 

(4.1)  ft  =  ux 

T£Th 

where  T  is  a  typical  (closed)  element  in  the  sinplicial  decomoosition  T.  ;  in 

h 

particular,  we  permit  nonsimplical  elements  near  the  boundary.  Let  denote  the 

linear  space  of  continuous  piecewise  linear  trial  functions  determined  by  T^: 

(4.2)  =  { x  «  C  (S)  :  X|Tis  linear  Vt  £  T^} 

Let  be  the  Ritz-Galerkin  H1  (ft)  projection  defined  by 

(4.3)  (E.V>,X)  i  =  Csf? » i  >  for  all  K  (  M.  , 

IT  (ft)  fC(ft) 

for  each  fixed  ^  e  H1  (fl) ;  here  we  use 

(4.4)  (*,<l>)  ,  =  <W>,Vi|»)  +  TT  I  *  /  1>  ' 

H  (ft)  L*2  («)  |fl|  SI  n 

which  defines  a  norm  equivalent  to  the  standard  H1  (0)  norm  in  the  usual  way  (cf.  [13]). 
One  sees  easily  that  preserves  (integral)  mean  values  over  ft  and  that  E^  is 

an  orthogonal  projection  in  H^(ft).  Finally,  define  the  finite  rank  approximation 
Th  of  T  by 

(4.5)  Th  =  EhoT  ,  Th  :  H_i(ft)  -*■  . 

Remark  4.1.  The  mapping  T^,  defined  by  (4.5),  is  not  the  natural  finite  rank 
projection  associated  with  T;  this  would  be  obtained  by  defining  E^  with  respect  to 

(4.6)  <4>,i|i)  =  (V*,VK») 

H^(ft)  u  (ft) 

leading  to  a  self-adjoint,  positive  definite  operator  on  It  is  therefore  some- 

2 

what  surprising  that  is  self-adjoint  and  non-negative  definite  on  L  (fl) . 

2 

I^mna  4.1.  T*h  is  self-adjoint  on  L  (Q)  and  satisfies 


(4.7) 


(T,g,g)  ,  >0,  for  all  g  i  L  (S2) 

h  l2  tn)  " 


Proof:  Since  T  is  bounded,  the  self-adjointness  will  follow  from  the  symmetr 

h 


of  T^.  This,  together  with  (4.7),  follow  from  the  following  identities.  Since  T 


is  a  Riesz  mapping, 


(T  g  ,g  )  ,  =  (VE  Tg  ,7Tg  ) 

h  1  2  L2  (Q)  h  1  2  L2  (tl) 


so  that,  by  (4.4),  the  mean  value  property  of  E^,  and  its  role  as  an  orthogonal 


projection  in  H  (£2) , 


(Thg,.gJ  ,  =  (E  Tg  ,Tg  )  ‘  -pT  /  /  Tg 

h  1  2  l2(£2)  h  1  2  Hha)  M  a  ^  1  a  2 


-  JnJ  v»2 


■  ('7EhT9l''7EhT92,L2(iJ) 


The  latter  quantity  is  non-negative  if  g^  =  g2  =  g  and  the  interchange  of  g^  and 


g2  yields  the  symmetry  of  T^ . 


Remark  4.2.  Our  basic  finite  element  hypotheses  are  the  following.  For  g  e  L  (Q) 


(i)  g  >  0  ■»  T.g  >  0  (Korovkin  Positivity)  ; 
—  h  — 


(4.8) 


(ii)  ||  (T  -  T  ) g  1 1  <  Ch  || g  || 

t  (n)  l  (a) 


holds  for  C  independent  of  h.  Note  that  (4.8ii)  is  a  routine  consequence  of 


(2.2)  for  m  =  2  since,  if  w  =  Tg,  then  =  T^g  satisfies  w^  =  E^w.  Ttie 


2  2 

order  h  approximation  properties  in  L  are  well-known  in  this  case  (cf.  [14]). 


We  are  now  ready  to  define  the  continuous  time  finite  element  approximation  Uh- 


The  approximation  is  based  upon  (2.13)  and  the  intrinsic  approximation  properties  of 


the  operators  T 


Definition  4.1.  The  finite  element,  continuous  time  approximation  :  1 0 , T^ ]  -*■ 


r  -  C  h  is  defined  to  be  the  unique  solution  of  the  initial  value  problem 
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(4.9i) 


!  "Vh't +  (uh  * r)  '  Th  r  -  V)(* '  V  i 


for  all  xer*CnMji,0<t<TQ  , 

(4.9ii)  Uh(0)  =  0  . 

Remark  4.3.  Of  course,  it  requires  an  argument  to  prove  that  (4.9)  possesses  a 
unique  solution.  Uniqueness  is  established  in  a  manner  similar  to  that  of  the  proof 
of  Theorem  2.4;  in  particular,  (cf.  (2.12))  it  can  be  proved  that 


(4.10) 


2  (ThlUl,h  '  Vh1,  Ul,h  "  U2,h'L2(n) 


♦  I 


l,h 


U  ||  ds  <  0 
2'  I,  (H) 


for  any  two  solutions  U,  .  and  U,  ,  of  (4.9),  where  (4.10)  holds  for  each 

1 1  n  2  $  n 

0  <  t  <  Tq.  Uniqueness  follows  from  (4.10)  and  (4.7).  An  existence  proof  can  be 

constructed  using  the  method  of  horizontal  lines,  wherein  the  time  derivative  in 

(4.9i)  is  replaced  by  a  backward  divided  difference,  quite  similar  to  the  technique 

in  section  three.  There  result  stationary  problems  which  are  the  gradient  formula- 

2 

tions  of  quadratic  minimization  problems  over  closed  convex  subsets  of  L  (2).  The 

step  function  (in  time)  sequence  indexed  by  M  ■  T/At  can  be  shown  to  be  bounded  in 
2 

L  (Q)  (actually,  the  finite  dimensional  affine  subspace  r  -  M^)  and  convergence 
of  a  subsequence  to  a  solution  of  (4.9)  may  be  established  by  standard  methods. 

Lemma  4.2.  Under  the  hypothesis  (4.8i),  the  relation  (4.9)  holds  for  all  x  £  r  -  C 
and  in  particular  for  x  =  u. 


Proof:  Since  Th  is  self-adjoint,  (4.9i)  is  equivalent  to 


(4.11) 


/  U(Uh,t  -  f](Thx  -  W  +  VUh  •  V(Thx  -  ThUh)}  >  0 


for  all  x  e  r  -  C  n  M^.  Me  rewrite  (4.11)  as 

(4.12)  /  {[(uh)t  -  f)  (;  -  Thuh)  +  vuh  •  vi;  -  Thuh)}  >_  0  , 

for  all  c  t  T^r  -  T^  (C  n  )  »  T^M^.  The  lemma  will  therefore  follow,  if  we  can 
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prove  that 


(4.13)  Th(r  -  C)  =  ThM^  , 
from  the  self-adjointness  of  T  .  Here  we  have  written 

(4.14)  =  (x  «  \  :  X  >  0}  . 


Since  T.  is  assumed  pointwise  non-negative  by  (4.8i),  (4.13)  will  follow  if,  for 
h 

v  e  L2 (ft) , 

(4.15)  4  €  1  T^}  =»  {4,  i  Thv)  . 

2 

Here,  the  orthogonality  in  (4.15)  is  understood  to  be  L  orthogonality.  To  verify 
2 

(4.15) ,  let  v  €  L  (ft),  let  e  and  suppose 

* 1  vs.  • 


By  the  proof  of  Lenrna  4.1  we  have 


0  =  (T  iji,  ip) 

h  L2(fl) 


l*V 


1  2 
L  (ft) 


so  that  T^iJ/  is  constant  on  0.  Thus,  using  the  aforementioned  proof  once  again  we 
have, 


(VT.v,VT.  0)  , 

h  h  L2  (ft) 


(T  v,  4/) 
h  L2  (ft) 


which  verifies  (4.15).  Hie  proof  is  completed. 

We  are  now  ready  for  the  major  result  of  this  section. 

Theorem  4.3.  The  finite  element  approximations  defined  by  (4.9)  are  convergent  to  the 

2  2 

solution  u  of  (1.10)  with  order  h  in  L  (D) ,  i.e., 

(4.16)  l|u  -  U  ||  ±Ch2 

1*  (D) 

for  some  constant  C. 

Proof:  Setting  z  =  in  (2.13ii)  and  x  =  u  ih  (4.9i)  (cf.  Lemma  4.2)  we 

have 
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(4.17) 


(i)  /  t(Tu)  +  (u  -  r)  -  TAr  -  Tf)  (u  -  U.  )  <  0 

r.  h  " 


(ii)  /  [<T  U  )  +  (U  -  r)  -  T  &r  -  T.  f]  (u  -  11)  >  0  . 

%  n  n  t  n  n  h  h  — 

Subtraction  of  (4.17ii)  from  (4.17i)  gives,  after  some  rearrangement, 

v2 


/  {Th(u  -  Vt  •  <«  -  V>  +  /  tu  -  V 


(4.18) 

Wri ting 

(4.19) 


ft 


/  (u  -  U  )  (T  -  T)u  +  /  (u  -  U.  )  (T  -  T.  )  (f  +  Ar) 

—  n  h  h  t  ‘  h  h 


I  (Th(u  -  0h)t  *  <u  -  0h,}  =  1  ^  /  Th(u  -  V  -  (u  -  0h)  , 

4*  44 


integrating  in  t  from  0  to  T  and  using  the  non-negative  definiteness  of  T 

0  h 

yield 

U.20)  llu  ~  U  I)2  £jl|u-u  ||2 

IT(D)  .  n  L2(D) 


+  F  IH*  -  V<Ut  -  f  -  AD  ||2 

L 2  (D) 


where  we  have  applied  the  standard  inequality  (3.10)  to  the  right  hand  side  of  (4.18). 
The  estimate  (4.16)  now  follows  immediately  from  (4.20)  and  (4.8ii)  applied  to 
g  ”  u  -  f  -  Ar.  The  proof  is  now  concluded. 
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vergence  rates.  The  following  explicit  results  are  obtained: 

( 3 i )  ||UM  -  u  ||  1  CAt  , 

L  (D) 

where  (/'  is  the  H1 (ft) -valued,  piecewise  linear  in  time,  interpolant  obtained 
from  the  horizontal  line,  fully  implicit  Euler  scheme  applied  to  (2)  with 
At  =  T  /M; 

(3ii)  ||U  -  u  ||  £  Ch2  , 

L  (D) 

where  U,  is  the  continuous  time,  finite  element  approximation  obtained  by 
h 

applying  an  integrated  version  of  (2)  to  a  translate  of  the  finite  dimensional 

trial  space  of  C°  piecewise  linear  elements.  The  approximation  scheme  used 

to  define  U,  appears  to  be  new. 
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